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ON THE CONTINUOUS RESONANT EQUATION FOR NLS 
I. DETERMINISTIC ANALYSIS 

PIERRE GERMAIN, ZAHER HANI, AND LAURENT THOMANN 


Abstract. We study the continuous resonant (CR) equation which was derived in [7] as the large-box 
limit of the cubic nonlinear Schrodinger equation in the small nonlinearity (or small data) regime. 
We first show that the system arises in another natural way, as it also corresponds to the resonant 
cubic Hermite-Schrodinger equation (NLS with harmonic trapping). We then establish that the basis 
of special Hermite functions is well suited to its analysis, and uncover more of the striking structure 
of the equation. We study in particular the dynamics on a few invariant subspaces; eigenspaces of the 
harmonic oscillator, of the rotation operator, and the Bargmann-Fock space. We focus on stationary 
waves and their stability. 


1. Introduction 


1 . 1 . Presentation of the equation. The purpose of this manuscript is to study the so-called con¬ 
tinuous resonant equation which was introduced by Faou-Germain-Hani [ 7 ] as the large-box limit of 
the cubic nonlinear Schrodinger equation in the small nonlinearity regime. This equation reads 

{idtu = T{u,u,u), {t,x)£RxR'^, 

[ u{ 0 ,x) = f(x), 

where the nonlinearity is defined by 


( 1 . 1 ) 


r(/i,/2,/3)(^)"= 


fi{x -L z)f2{Xx'^ -L z)f^{x + Ax-*- -L z) dxdX 


/i(x“^ -|- z)f2{Xx 2;)/3(x-*- -|- Ax -|- 2:) dx dX, 


for any z G (if x = (xi,X2) we set x-*- = (—X2,xi)). While the above formula seems mysterious at 
this stage, it can be thought of as an integration over all rectangles for which 2; is a vertex. Indeed, 
the points 2;, x -|- 2;, Ax-*- -|- 2;, x -|- Ax-*- -|- 2; form a rectangle in and this yields a parameterization 
of all rectangles which count 2; as a vertejiQ- 
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This expression can also be reformulated using the unitary group as was observed in [7]. 

Note that the definition of T above is slightly different but equivalent to that in [7| as we explain 
in Section ri.dl Here we will show that T can also be reformulated using the semigroup ) 

(see Lemma [22] below). The key ingredient is the so-called lens transform which is an explicit for¬ 
mula (given in (12.5M which links to ). This suggests that the harmonic oscillator 

H = — A -|- |xp will play a central role in the study of (ICRh . 


Defining 


def 


^(/l,/2,/3,/4) = (r(/l,/2,/3),/4)L2 

= f [ [ fi{x + z)f2{Xx^ + z)f3{x + Xx-^ + z) f4{z)dzdxd\, 

Jr ./r2 Jr2 

it is easy to check that the (ICRh equation derives from the Hamiltonian 

^(/)"= SifJJJ) 

given the symplectic form oj{f,g) = —4Jm(/ , L^(M^) (this follows easily from the symme¬ 

tries of T). In other words, (ICRD can also be written 

tj 2 df 

Important quantities conserved by the flow of the above equation (we shall come back to them) are 
the mass M and angular momentum P: 


f 

Jr^ 


m"!/ / InP 


and P 


/ i{x X V)uu, 

./r2 


where x x V = X 2 dxi — xid. 


■'X2 ■ 


1.2. Physical and mathematical relevance. The (jCRp equation has rich dynamics and can be 
studied in its own right, but it also plays a role in the description of the dynamics of the usual cubic 
NTS - with or without potential - in various situations, which we summarize here: 

• It was derived in [7] as a weakly nonlinear, big box limit of the cubic NLS 

idtu — Au = |npn 

(here, we consider the focusing case, but the defocusing case leads to the same picture) on 
a periodic box of size L; equivalently, it appears as the limiting equation for high frequency 
envelopes of solutions of NLS on the unit torus T^. To be more specific, setting the above equa¬ 
tion on the 2-dimensional torus of size L, and prescribing data of size e, it is well-approximated 
by (CR) on very long time scales (much longer than 1? j^'). 

• We will prove in the present paper that (ICRp can also be derived as a small data approximation 
of the 2-dimensional cubic NLS equation with harmonic trapping, a.k.a. Hermite-Schrodinger 
equation 

idtu — Au P |xptt = \u\^u 

This model is widely used in several areas of physics from nonlinear optics to Bose-Einstein 
condensates [16], but its relation to the dynamics of (jCRp seems to be new. To be more 


( 1 . 2 ) 
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specific: set H = —A + |xp, n„ the projector on the n-th eigenspace of H, and f = e 
Keeping only the totally resonant part of the nonlinearity in (ll.2p gives the equation 

idtf= Iln,[{UnJl){Unj2)(n^3)). 

ni,n2,n3,n4>0 

n-i+n2=n3+n4 

We will prove that the above right-hand side is, up to a multiplicative constant, equal to 
T(/, /, /). Thus the totally resonant part of NTS with harmonic trapping is identical to (ICRh . 
This implies that (ICRh approximates the dynamics of (II.2D for large times in the small data 
regime, as we illustrate in Theorem 13.11 

• The equation (ICRD appears as a modified scattering limit of the cubic NTS on with harmonic 
trapping in two directions. Therefore any information on the asymptotic dynamics of (ICRD 
directly gives the corresponding behaviour for NTS. We refer to Hani-Thomann [H] for more 
details and concrete applications. 

• When restricted to the Bargmann-Fock space, which is given by functions which can be 

-Ml 

written as the product of a holomorphic function with the Gaussian e 2 ^ the equation (ICRp 
coincides with the model known as Lowest Landau Level [iiiaDS], used in the description of 
rotating Bose-Einstein condensates. 

The (jCRp equation can be considered as a model of NLS-like equation without dispersion. A lot 
of attention has recently been paid to such equations, at least starting with the work of Colliander, 
Keel, Stafhlani, Takaoka, and Tao [5] and subsequent works on the growth of Sobolev norms [121 
13]. There, the dynamics of resonant systems like (CR) arise either as approximating or asymptotic 
dynamics for the original NLS model (see also M)- Another important instance of zero-dispersion 
Hamiltonian equations comes from the work of Gerard and Grellier mm who studied the so-called 
cubic Szego equation on (this latter equation is obtained by replacing H^ in T (see Lemma [2j3l) by 
the projection on for x G §^). We refer to Pocovnicu [nncn] for the study of Szego on M. 

Despite the absence of a linear part, and of the corresponding dispersive effects, dCR]) is well- 
posed in L^(M^) (see m) which is remarkable for a zero-dispersion equation with a zero-order trilinear 
nonlinearity (in comparison, the Szego equation is well-posed in for s > 1/2 and the result is 
sharp). In fact, the nonlinearity has a hidden smoothing property coming from Strichartz estimates 
of the original NLS model, and this compensates the lack of linear dispersion (see Lemma l2.2p . In 
this direction, let us also mention the dispersion managed Schrodinger equation, which is obtained by 
averaging a nonlinear Schrodinger equation with varying dispersion. This latter equation has a similar 
structure as dcEl). We refer to [Sam] and references therein for more details. 

In [TU] we undertake the study of (IGRD with random initial conditions. We exhibit global rough 
dynamics (for initial data less regular than L^), and we construct Gibbs measures which are invariant 
by this flow. This is related to weak turbulence theory, which is often considered to occur in the 
statistical regime where the phases of the Fourier coefficients are initially uncorrelated. 

1.3. Known results. Here we recall some important properties of the (IGR.D and the operator T, 
which were proved in [7]. We first clarify the slightly different formulation of the operator T we use 
here and that derived in [7]. In addition to reversing the roles of /2 and /s, the main difference between 
the two definitions is that the A integral in [7] is over [—1,1] instead of M. Denoting by Tvgf-ip] the 
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same operator as T with the A integral taken over [—1,1] instead, one can easily show by a change of 
variable that: 

rA6[-i,i](/i,/i,/3) = ^r(/i,/i,/3). 

In particular, one can use either formulation to define the (jCRI) equation in which fi = f 2 = fs- 

Now we recall that the operator T is bounded from L'^ x L'^ x to and also from ^ ^ 

to where L°°’^ is given by the norm ||/||^oo,i = |||x|/||loo. 

This implies immediately that (jCRD is locally well-posed for data in or Using the conser¬ 

vation of the norm (the mass M), one obtains global well-posedness for data in (and afterwards 
in any Sobolev and weighted space). 

Gaussians play the role of a ground state for the equation, since they minimize the Hamiltonian £ 
for fixed M. This variational characterization leads to orbital stability in L?' (up to the symmetry 
group of the equation, which will be recalled in Section EH); it also holds in n (where is 
given by the norm ||/||£ 2 .i = ||(a;)/|| 2 , 2 ) by a different argument. 

1.4. Obtained results. 

1.4.1. The basis of special Hermite functions. Recall that there exists a Hilbertian basis of L^(]R^) 
known as the special Hermite functions {(pn,m}, where n € N, m £ {—n, 2—n,..., n —2, n} which diago¬ 
nalizes jointly the harmonic oscillator H = — A-|-|xp and the angular momentum operator L = ix x V: 

H^n,m — 2(71 -|- T)'pn,mi LTn,m — '^Tn,m 

(see Section[5]for a more detailed presentation). We show in Section [5] that this basis is very well-suited 
to decomposing the trilinear operator T: this follows from the formula (proved in Proposition 15.2h 

T (v^m ,mi) Tn2,m21 ,m3 ) — ^(v^ni ,mi ) Tn2,m2 ) ) ¥^ 714,7714 )¥^n4,m4 j 

where 714 = ni -|- n 2 — 713 and 7774 = mi + m 2 — 7773 . This formula has the immediate consequence that 
all the special Hermite functions are stationary waves. Furthermore, it also implies the dynamical 
invariance of all the subspaces of the form 

Spanj ;^2 {ipn,m, such that an + Pm = 7 mod <5} or Spanj ;^2 {ipn,m, such that an + Pm = 7} 
where a, P, 7 and 5 are natural numbers. 

1.4.2. Invariant subspaces. We investigate further the dynamics on particularly relevant or natural 
examples of the invariant subspaces which were just described: 

• Sections 16.1116.41 are dedicated to the analysis of the dynamics on the eigenspaces of the har¬ 
monic oscillator H: for some ttq G N, Span £2 {ipno,Tn, where m G {— 77 o, 2 — no, ..., 779 — 2, 779 }}. 

• Section [7] focuses on the eigenspaces of the rotation operator L: for some 7779 G Z, 

Span ^2 {yjn,mo, where n G {jmgl, \mQ\ + 2, 17779 I -h 4,... }}. 

• Finally, in Section [HI we study the equation on the Bargmann-Fock space which can be seen 

as Span 2,2 where n G N}. 






5 


1.4.3. Stationary waves. The most classical notion of a stationary wave is given by a solution of the 
type where is a fixed function, and cu G M. More elaborate waves are of the type 

where Rq is the rotation operator in space around 0 of angle 9, and a, to real numbers. For reasons 
that will become clear, we call the former M-stationary waves (or simply stationary waves), and the 
latter M + aP-stationary waves. 

In the invariant subspaces detailed above, we give examples of stationary solutions, and try to 
investigate their stability. Understanding the full picture - even finding all stationary solutions - 
seems a daunting task, but we obtain first results in this direction. 

Finally, in Section [H we prove general theorems on decay and regularity of stationary waves of (ICRD . 
More precisely, we show that any M-stationary wave in is analytic and exponentially decreasing, 
while, roughly speaking, M + aP stationary waves belong to the Schwartz class as soon as they are a 
little more localized, and smooth, than L^. 


1.5. Notations. We set 

• For X G x"*" is the rotation of x by ^ around the origin. 

• For X G M^, (x) = yl + lxp; similarly if x G M. 

• (/, 5 ) 1 , 2 (r 2 ) 42 fg. 

• {0 = 1 ( 0 '= ^ [ /(x)e-“«dx. 

• N is the set of all non-negative integers (including 0). 

• H '^= — A -|- |xp is the quantum harmonic oscillator on M^. 

• L i{x 2 dxi — Xidx 2 ) is the angular momentum operator. 

• is the Sobolev space given by the norm 11/14'’*=^ II(?)^/IIl 2 . 

is the weighted L^-space given by the norm ||/||l 2 ,s '^= ||(x)^/ 42 . 

• P* = P* n is the weighted Sobolev space given by the norm ||/||h'> II( 2 ^)*/IIl 2 -|- 

II(0^/IIl 2. It is classical (see (|3.ip l that = {u G L^, s.t. G 

• Re is the counter-clockwise rotation of angle 9 around the origin. 

• Sxu = Att(A-) is the scaling. 

• is the orthogonal projection on the eigenspace = {u G L^(M^), Hu = 2(n -|- l)u }. 

In this paper c, C > 0 denote universal constants the value of which may change from line to line. 
For two quantities A and B, we denote A < B A < CB, and A ^ B A < B and A> B. 


2. Properties and symmetries of T and £ 

2.1. Various formulations of P and £. As observed in [7] we have 

Lemma 2.1. The quantities P and £ are invariant by Fourier transform: 

-^(r(/l,/2,/3)) =r(/l,/2,/3) and T(/i,/2,/3,/4)=^(/i,/2,/3,/4). 


( 2 . 1 ) 
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Proof. Using Fourier inversion and the identity f dx = gives 

^T{f,g,h){^) = -^ f f [ e~^^^f{x + z)g{Xx'^ + z)h{x + \x^ + z)dxdXdz 
2vr J]R2 Jr Jr2 


Ihvr^ 



gi^(-^+a+/ 3 - 7 )gix(a-A /3 +A 7 -/) f(a)g{P)h{j) da d^ dj dx dX dz 

dm \p^) dp 


Changing variables to 77 = ^(^ — /3)-*- gives 

J^Tif,g,h){^) = [ [ f{C + g)g{^ +Xg-^)h{g +Xg-^+^)dgdX, 


which is the desired result. 


□ 


The next result shows that £ can be related to the Strichartz norm associated to the linear 

flows e®*®^ and e“®*^. 

Lemma 2.2. The following formulations for £ hold 

(2.2) £{fl,f 2 ,h,U) = 271 / / {£^^fim^f 2 )(m 3 )(^ 4 )dxdt 

Jm. JR 2 

TT 

(2.3) " . 


= 271 / / (e-®^"/i)(e-®‘"/2)(e-®i^^/3)(e-®i^^/4)dTdt. 

J-f Jr2 


Therefore we have 


dt 


(2.4) 


r(/i,/2,/3) = 271 / e-®*^[(e®*^/i)(e®*^/ 2 )(e®‘A/ 3 ) 

JIR L 

= 277 |%®‘^[(e-®‘^/i)(e-®‘^/2)(i^^) 


dt. 


Proof. Using Fourier inversion and the identity ^ fj ^2 dx = 6^=0 gives 
mflW ^^ f 2 )(^ 3 )(^ 4 )dxdt 




1 


(277)4 

1 


g-*i(|a| +I/II -|7| -|5| )e^^(<^+f^-^-^)g{a)f2{P)Ml)fAi5)dadpdgd6dxdt 


+I/3I -|«+/3-<5| )f,(a)MP)Ma + P-d)h{d)dadpd5dt. 
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Changing variables to S = z,a = z + x ,/3 = z + Ax-*- + //x and resorting to the identity ^ dy = 

6^=0 yields 


A = 


1 



(277)2 
1 

= f2, fs, Ia) = f2, fs, Ia) 


+ x)f2{z + Ax-*- + yLx)f2,{z + x + Ax-*- + yx) /4(z)|xp dx dzdXdydt 



/i(x + z)f2{z + Ax-*-)/3(z + X + Ax-*-) fAiz) dx dz dX 


2tt 


which gives (|2.2I) . Let us now prove (I2.3p . Let / G L2 (m 2), and denote v{t,-) = e f and 
u{t, ■) = e®*'*^/. Then the lens transform gives (see for instance [22]) 


(2.5) 


u{t, x) = 


Vl +4t2' 


arctan(2t) x 

2 ’ vTW ■ 


3/ QirctOiii ^ ^ 

We first make the change of variables y = =, then r =-. This gives 

V1 + 4t2 2 

£{fij2,h,fA) = I 

= 2 tt [ [ {e~"^f2)ie-^'^^f3){e-^^^U)dydr. 

J-IL Jr2 


The relations for T are obtained using that (T(/i, /2, /s), /4 )l2(r2) = £{fi,f2, /s, /4)- □ 

We are now able to prove the following result 
Lemma 2.3. The following formulations for £ hold 

£{flj2,f3,fA)=7r^ [ ^{Iln^fl)innj2){n^3){K:jA)dx. 

ni+n2=n3+n4 


Therefore we have 

r(/i,/2,/3) = 7r2 n,,((n,ji)(n,,/2)(n;:;:^)). 

ni+n2=n3+n4 


Proof. We compute £ in (12.31) for the eigenfunctions of H. Therefore we assume that Ilnjfj = fj and 
then 

^(/l,/2,/3,/4) =277 e-2®(®®l+®®2-®®3-4)irfi f f.fjjfdx. 

J-f JR2 

But now we use that n„^/(—x) = (—l)®®7n„^/(x), thus J^2 fif2f3fAdx = 0 unless rei + n2 — n3 — 77.4 
is even, which in turn implies Q-'^dm+n2-n3-n4)tj^^ _ + n2 — — 774). □ 
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2 . 2 . Symmetries of T and conservation laws for (ICRD . 

Lemma 2.4. The commutation relation 

( 2 . 6 ) Q(r(/i, / 2 , /a)) = r(Q/i, / 2 , /s) + r(/i, Q/ 2 , /a) - r(/i, / 2 , g/a) 

/loWs /or t/ie (self-adjoint) operators 

g = 1 , X , |xp , iV , A , H , L = ix X V , i{x • V + 1) 
for all /i,/ 2,/3 G T^{Q), where T>{Q) denotes the domain of Q. 

Proof. First observe that it suffices to prove the commutation relation for /i, / 2 , /a sufficiently smooth, 
and then argue by density and boundedness of T. For g = 1, x, |xp, this follows easily from the 
definition dEU) of T, in particular the fact that the arguments of fi, f 2 , /a satisfy 

z + {x + \x^ + z) = {x + z) + (Ax"*" + z) 

\z\^ + |x + Ax"*" + = \x + z\ + lAx"*" + 

Using m and arguing similarly gives (ESI) for P = zV and A. Combining |xp and A gives (|2.6p for 
Q = H. Finally, to obtain this commutation relation for g = zx x V and z(x • V + 1), define R\ to be 
the rotation of angle A around the origin, and Sx the scaling transformation Sxu = Azz(A-), observe 
that 

RxT{u, u,u) = T{Rxu, Rxu, Rxu) 

SxT{u,u,u) = T{Sxu,Sxu,Sxu) 

and differentiate in A. □ 

Corollary 2.5. If Q is as in Lemma \2.4[ and /i,/ 2,/3 £ L^(M^) are such that Qfj = Xjfj, with 
Xj G C and j = 1, 2,3, then 

g(r(/i, / 2 , /a)) = (Ai + A 2 - A 3 )r(/i, / 2 , /a) 

and for s G M, 

e“®r(/i,/2,/3)=r(e“'^/i,e“0/2,e'^«/3) and f) = £{f). 

In particular, if /i, / 2 , and /s are eigenfunctions of Q as in Lemma 12^ then so is T(/i, / 2 , /a)- This 
corollary hints towards examining 7” in a basis which simultaneously diagonalizes both H and the 
angular momentum operator L = i(x x V) = i{x 2 dxi — this will be done in the next section. 

Lemma 2.6. If Q is an operator so that for all f G 5(M^) 

Q{r{f, /, /)) = 2 T{Qf, /, /) - T{f, /, g*/), 

then 



is a conservation law for dcED. Applying this to 

g = 1 , X , |xp , iV , A , H , L , z(x • V + 1) 
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gives the conserved quantities 
M = I \u\\ I , / IxiVP , / ^VuW , / I V„|^ . I Huu.P = 


= / Luu 


(which are real-valued since we chose Q self-adjoint) 
Proof. We compute 

(L 

dt 


(2.7) 


f {iQu)u = {QT{u),u) - {Qu,T{u)). 
Jr'2 


By assumption and the symmetries of £ 

{QT{u),u) = 2{T{Qu,u,u),u) — {T{u,u,Q*u),u) 
= 2{Qu,T{u))-{T{u),Q*u), 

which implies that {T{u),Q*u) = {Qu,T{u)), and yields the result by (12.71) . 


J i{x ■ V + l)uu 


□ 


The relation which has just been established between operators commuting with T, symmetries 
of T and T, and conserved quantities of (ICR| is of course an instance of the Noether theorem. We 
recapitulate below the obtained results (with A G M). 


operator Q 
commuting with T 

conserved quantity 
f Quu 

corresponding 
symmetry u i-A e^^^u 

1 

j>l^ 

U eA e^^u 

Xl 

f Xl\u\'‘^ 

u i-A e^^^^u 

X2 

f X2\u\'^ 

u i-A e^^^^u 


/|xp|«|2 

u i-A 


J idx^uu 

ue^ u{- + Aei) 


J idx^uu 

ut-^ u{- + Xe 2 ) 

A 

TWW 

u i-A e^^^u 

H 

f Hull 

u i-A e'^^^ u 

L = i{x X V) 

J Luu 

u 1 —>■ u{R\x) 

z(x • V + 1) 

J i(x ■ V + l)uu 

u 1 —>• Au(Ax) 


3. Approximation of NLS with harmonic trapping by (CR) 

A first consequence of the findings of the previous section is the following theorem, which states 
that (jCR.j) approximates the dynamics of (|1.2p in the small data regime. For s > 0, we define the 
Sobolev space based on the harmonic oscillator "H® = jn G L?, s.t. G endowed with the 

natural norm HrH-rs = \\H^/‘^u\\i 2 . By [241 Lemma 2.4], we have the following equivalence of norms 


(3.1) 
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Theorem 3.1. Let s > 1 and suppose that u{t) is a solution of (|1.2I) and f{t) a solution of (ICRI) 
with the same initial data uq. Assume that the following bound holds over an interval of time [0,T] 

\\f{t)\\w<B for all 0<t<T. 

Then there exists a constant C > 1 such that 

\\u{t)-e^*^f{^)\\ns(M. 2 ) <C{B^ + for all 0<t<T. 

In particular, if B is sufficiently small and 0 < t <C B~^logB~^, then 

ll«(*) - e‘‘"/(^)ll«-(K») < 

We remark that B can be made small by taking the initial condition sufficiently small in 
Proof. We start with some notation: we write oj = ni + n 2 — — n 4 and 

r'(/i,/2,/3)= n,,((n„ji)(n„j 2 )(n;;;;^)), 

ni,n2,?T-3,n4>0 

M(/i,/ 2,/3) = e**“n„,((n„Ji)(n„,/2)(rra)) = 

ni,n2,n3,n4>0 

mij2j3)= Y1 e**-n„,((n„ji)(n„j2)(Tw^)) =a4(/i,/ 2,/3)-r'(/i,/2,/3). 

ni,Ti2,n3,n4>0 

uj^O 

Note that T' only differs from T by a constant multiplicative factor (see Lemma l2..Sp . Recall that 
for s > 1, we have C and that is an algebra. Then by (|3.ip . it is easy to 

check that is also an algebra when s > 1. Therefore, using the boundedness of the operator 

giiiL g_i^s (uniformly in t), we obtain the boundedness of the trilinear operators T' and V from 
X LL^ X LL^ to LL^ for any s > 1. 

Let g{t) = e~'-^^u{t) and f{t) = f{^). The equation satisfied by g{t) and / are the following: 

idtg{t) = Nt{g{t),g{t),g{f))-, idj = T {1,1,7)] g{0) = f{0) = uq. 

Now we write the equation for g{t) as follows: 

idtg =T'{g,g,g) + ^ (^(n„i5r)(n„25)(n„35)^ 

ni,n2,n3,n4>0 

(3-2) =T'{g,g,g) + dt ^ ——— Un^\^{Un^g){Un 2 g){'nn 3 g)j 

ni,n2,?T-3,n4>0 

uJ^O 

_^ / _ s 

(3-3) — 'y ^ Iln4^dty{Llni9){B-n2 9){Bn3 9)j 

ni,n2,n3,n4>0 

CxJf^O 

=T'{ 9 , 9 , 9 ) +dtA + D, 
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where A is the sum in (13.2p and D is given in (13.3p . Now denote by [t] the largest integer smaller 
than t and notice that 


E 


gito; _ I 


ni,n2,n3,n4,>0 
rt 


lUl 


nn4 ( (nni9)(nn,25)(nr),35) ) ~ 


ni,n2,n3,n4>0 

uj^O 


e^^^UnA{Un,g)iUn,g){Un,g)] dr 


[ 'PT{g{t),g{t),g{t))dT. 


This allows to estimate 


II^IIh- A\\dtg\\w\\g\\w ^ llffllw 


Now let e{t) = g(t) — f{t), and assume as a bootstrap hypothesis that ||e(t)||'^s < B. Then the 
equation satished by e{t) can be written as 

dte = T'{e, g, g) + T'{J, e, g) + T'(f, f, e) + dtA + D. 

This gives that 

\\eit)\\n^ < CiB^ [ ||e(s)||H»ds + C^B^ + tB^- 
Jo 

The result now follows by Gronwall’s inequality. □ 


4. Regularity and decay of stationary waves 
Recall the following conservation laws 

M{u) = I |ttp, P{u) = I Luu. 

As we mentioned in the introduction, M + aP-stationary waves read where a and oj 

are real numbers, and 99 is a fixed function. In particular, M-stationary waves are simply of the type 
Notice that, in degenerate cases (if is an eigenfunction of L), a given solution can be an 
M + aP-stationary wave for al0 a. Variationally, M + aP-stationary waves can be characterized 
as critical points of £ under the constraint that M + aP takes a fixed value. The Euler-Lagrange 
equation reads 

(4.1) u}(p + au } L(p = T {^ p , i ^,( p ). 

In the following subsections, we prove polynomial or exponential decay, in the physical space or 
Fourier variable, for large classes of such stationary waves, requiring that they belong to L^, or slightly 
more. 

On the one hand, the condition of belonging to is sharp in order to obtain decay or regularity. 
This can be seen through the example of j^, which is a stationary wave in the weak-L^ space 
but not in any Sobolev or weighted space of positive index. On the other hand, all the examples we 

2lt would be natural to define P-stationary waves as waves of the type R-auit^P- However, we could not find an 
example of such a wave which is not at the same time an M-stationary wave. We conjecture any P-stationary wave is 
also an M-stationary wave. 
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know of M + aP stationary waves are products of Gaussians and polynomials, suggesting a possible 
improvement of our results. 

4.1. Polynomial decay. 

Theorem 4.1. Assume that a G M, with a = Q or a~^ ^ Z, and let be an M + aP stationary wave. 
Then 

(i) If ip belongs to for some e > t), it belongs to for any s > 0. 

(a) If ip belongs to , for some e > t), it belongs to for any s > 0. 

(in) If p belongs to n for some e > t), it belongs to the Schwartz class. 

Proof. By invariance of the Euler-Lagrange equation (14.1|) under the Fourier transform, it suffices to 
prove (i). Using the fact that (Id+aL)“^ is bounded on - since L commutes with radial weights 
- it will follow from the repeated application of the following proposition to (|4.1I) . □ 

Proposition 4.2. For a > 0, and 5 < the operator T is bounded from to . 

Remark 4.3. This proposition also implies a very strong smoothing effect for the equation under 
study. 

Proof. By duality, it suffices to prove that, for / and g in with norm 1, which are fixed from now 
on, and which we assume to be non negative, 

(r((x)-/, {x)-^f, {x)-^f),{xr+^g) < 1, 

or in other words 

[ [ [ K{x, z, X)f{z + x)f{z + Xx'^)f{z + X + Xx'^)g{z) dx dz dX < 1, 

J-l ilR2 JlR2 

where 

K{x, z, X) ^ x)^{z + Xx-^)^{z + X + Ax-*-)'^ ’ 

Step 1: bounds on K. Localizing dyadically z and z+x+Ax"*-, let us assume from now on that {z) k, 2^ 
and (z + X + Ax-*-) ~ 2^, where j and k are integers. 

First, it is clear that K <1 unless \z\ is large. The identity 

+ \z P X + Ax"*"!^ = \z P xf P\z P Ax"*"!^ 

implies that min(|z + x], \z + Ax-*-|) > \z\. If in addition \z P x P Ax-*-| > \z\, the assumption 6 < a 
entails that K < 1. 

Thus, K P 1 unless 2^ is large and 2^ < So2^, for sq > 0 chosen sufficiently small. Let us assume 
for now that these two conditions hold. 

There holds then \z + Ax-*-| Ri \z\. Indeed, \z P Ax-*-| ^ \z\ would imply that |Ax| « l^j, and then 
\z P X P Ax-*-! > |x| — 1^; + Ax-*-! > \z\ (since |A| < 1), which contradicts 2^ < eo2P 

Next, if lAx-*-! > 2eo2^, then |z + x| > lAx-*-! — \zpxp Ax-*-! > £o2P Combined with \zp Ax-*-! « \z\, 
and, once again 5 < a, this implies K <1. 

Therefore, K <1 unless 2^ large, 2^ < £o2^, and |A||x| < 2eo2L Under these three conditions, 
\z\^2\ |x| Ri 2-^', \zPXx^\^2^, [z P X P Xx^) ^ 2^ < £q2P 
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One can then estimate 

Iz + xl > |A||3:| — \z + X + \x^\ > |A|2'^ provided |A|2'^ > Co2^, for a constant Cq. 

When the previous estimate is valid, one finds by a straightforward computation 

2j{cr+S) 

\K\ < -< 

and this last bound is 0(1) if |A| > 2“(^+02'?a. 

Summarizing: K <1 unless 2^ large, 2^ < eo2^, and |A| < max(C'o2^“-^, 2“*^^'''02-? a). If on the other 
hand these three conditions are satisfied, we use the trivial bound 


(4.2) 


Step 2: Decomposition of T- Define 


2J0+<5) 

\K\ < . ^ < 2^^2-’^^. 


Xj,k{x,Z,X) l2><{2)<2>+ll2'=<{z+a:+Aa:-L}<2'=+ll^^^3^(^^2fc-^2-(fc+j)2ji)’ 
X{x,Z,X)= ^ Xj-fc(x, 2 :, A), 


and decompose 


2'“<£o2J 


K{x, z, X)f{z + x)f{z + Xx'^)f{z + X + Xx'^)g{z) dx dz dX 


+ jjj ...(i-x{.,,.X))d.dzdx^Ji + n. 

Using that |ifr| < 1 on the support of the integrand of II, we obtain immediately, by boundedness 
of T from (U^)^ to L^, 

//<|(r(/,/,/),<7)| <1. 

As for I, the bound (|4.2I) gives 


/< 2^V^'^|(7^-fc(/,/,/),5)|, 

2'=<eo2:'' 


where 


Tj,k{f,g,h){z) = 


1-1 ■ 


f{z + x)g{z + Xx'^)h{z + x + Xx'^)xj,kix, z, A) dx dX. 


Bym Proposition 7.7], the operator Tj^k ■ (A^)^ has an operator norm such that \\Tj^k\\{L^)3^L^ < 

max(C'o2^“'^, 2“(^+'^)2'^A). This implies that 

I< 2'?'‘52-^'"max(Co2^-^2-(^+^)2^l) < 1 

2'= <£o2J 


(since 6 < ^^), which concludes the proof. 


□ 
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4.2. Exponential decay. The purpose of this section is to prove analyticity and exponential decay 
of M—stationary solutions of (ICRI) . Indeed, we will show the following: 

Theorem 4.4. Any solution of the equation 

(4.3) uup = 

satisfies 


for some fi > 0. 

As a corollary, one gets that ip is analytic in a complex strip {( 21 , 22 ) € : |/m 2 i|, |/m 22 | < 

/i/2}. We follow the elegant proof of Erdogan, Hundertmark, and Lee [6] (see also|llj'l where the 
corresponding result is proved for dispersion managed solitons in ID, with the difference that cubic 
NLS in 2D (from which (jCE.p is derived) is critical as opposed to subcritical. This is reflected 

in the key estimate (j4.5p where the “dimensionless” gain of (^) is much weaker than the 
decay in the ID case (both gains are dictated by scaling). 


We start with some notation: Let 


Ff,,eix) = p 


\x\ 


1 + £\x\ 


and define the following weighted versions of the operators £{fi, f 2 , f 3 , fi)- 
(4.4) h, f3, fi) = £{e~^'^’^fi, e~^^’^f2, h, 

We shall need the following key estimates on the multilinear operator 
Proposition 4.5. Suppose that /i ,...,/4 G L^(M^), then the following estimates hold uniformly v 


p,£ 


(^) I^M..(/l,/2,/3,/4)|<nil/^lli^- 

i=i 

( a ) Suppose that for some i,k G {1,2,3,4}, one has supp fe C B{ 0 , M) and supp fk C B{ 0 , Nf 
M <C N, then 

//Vf\V2 4 

(4.5) i^M,e(/i,/2,/3,/4)i < f 


with 


Proof. The proof of part (i) will follow from the bound 

(4.6) - + A 2 ) - F^,,(e + A 2 + 2 ^) - F^,,(e + 2 ^) + F^,,(0 < 0 

for any p,£. Notice the elementary inequality: for all e > 0, r/i, r /2 in 

^ \m + m\ ^ \vi\ + ^ 2 ! ^ l??il 1 ^ 2 ! 

1+ e|T/i+r/2| “ 1+ e(|?/i| + |?/2|) “ 1+ e|i/i| 'i-+ £\m\' 

We get (14.6p by applying (14.711 twice, since ^ ^ + A 2 — (^ + A 2 + z^) + (C + z^)- 
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The proof of part (n) follows by using a refined bilinear Strichartz estimate: Indeed, without loss 
of generality, we can assume that fi > 0, and then 

/2, /s, /4)| < £{fl, /2, fs, h) 


<2,r||(e‘“/i)(e“/2)(e‘“/3)(e“/4)||t._ 



where we used in the third inequality that ||e**'^/e**^ 5||2 = ||e**^/e**^ 5 || 2 , in the next to last step 
the bilinear refinement to the Strichartz estimate in 2D (cf. Bourgain [3]) and in the last 

inequality the standard ^ Strichartz estimate. □ 


Remark 4.6. Since T(/i, / 2 , /s, fi) = £{fi, f 2 , fs, fi), the same estimates above for hold for the 
operator 


£fi,e{h,f2, /s, fi) 


£,,,eifl,f2,f3ji) 


where P = —iV, and with the corresponding support assumptions in part (ii) of Proposition 
imposed on fi and fk- 


Proof of Theorem \4.4\ We start with some notation. Assume for simplicity that ||y?||L 2 = 1. Let 
M > 2 to be fixed later, and denote 

f<{x) ■■= fix)X\x\<M, Uix) ■■= fix)XM<\x\<M2, f>{x) ■= f{x)X\x\>M2- 

We would like to find M,p, > 0 such that < 1 uniformly in e > 0. Since ||e^''’®/<||L 2 , 

fr^Wi^ < II/IIl 2 , the hardest part is to bound 11^2, which is the aim of the three 

steps below. 


Step 1: Let ijj := e^>^'^tp. We first obtain an estimate on ||V'>||l 2 . We start by multiplying (14.3p by 
e^^f^’^^^'lx\x\>M'^'i^, gives after integration 




£^l,e y?, (/?>). 


Therefore 


IIl 2 = i/>, i}>) 

= £ii,e{f’< + V'- + V’>,'0< + V’~ +'0>,V’< +"0- + 
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Expanding the terms in this sum, we bound terms that are cubic and quartic in using Propo¬ 
sition [32] by ||V'>lli 2 and ||i /'>||^2 respectively. Terms that involve at least one copy of are 
bounded using part (n) of Proposition 14.51 This yields 

w||i/’>iii2 < 

^ llV'>lli2 + IIV’>lli2 + {m ^'^^||i/^<||| 2 + ||V’~lli2^ Il'0>lli2 + ^'^^||i/^<|||2 + ||V’~lli2^ IIV’>IIl2 

< Il^>lli2 +e'^^^||V'>||i2 (m-V2 + ||^^||2^ +e3MM^ ('^- 1/2 ^ ||^^||3^^ ||^^||^^_ 

Dividing both sides by ||'0 >||l 2 we obtain 
(4.8) w||'0>||l2 < 

c(||V’>|li2 +e'^^^||iA>||i2 (m-V2 + ||<^^||2^) (^M-V2 + ||^^|||,)) . 


Step 2: Let u := ||'0>||l 2 and choose fi = M Estimate (14.81) above translates into 


(a; - CoM- 1/2 - Co||c/p..||i 2 ) u - Cqi/' - < Cq (m-^s + . 


for a new constant Cq > 0 

Let G(p) = — Cqu^ and denote by Umax the point at which u 1 —)• G(u) achieves its 

maximum on [0, 00 ). Let uq = Vmaxl‘^ and choose M large enough so that the following two conditions 
are satisfied: 


(1) C'o(M 1/2 + ||y?.^||^ 2 ) < min(a;/2,G(po)), 

(2) II</2~||l2 -I- ||v2>||l2 < 

This is possible by the monotone convergence theorem. 
Therefore, (14.8p gives 


G(||e^-^:/.>||^2)<G(po) 


for all e > 0. 


Step 3: Since 11^2 is continuous in e, we obtain that : e > 0} is contained in only 

one of the two connected components of G“^(— 00 , G(po))- Now 

l|e'^''’V>llL2 < ||e^i^||Loo||(^>||^2 < 2||y?>||L2 < uq, 

(using our choice /r = M“^ < 1/4 above), which gives that ||j ;,2 : e > 0} is contained in the 

component [0,uo], which means that 

||g^’M.e(^> 11^2 < pQ) for all e > 0. 


We can now complete the proof. Taking the limit e —0, we obtain by monotone convergence that 
||eA‘|a;|y ,>||^2 < uq, which implies in turn ||e^l^l(/9||j;^2 < 00 . A parallel argument using Remark 14.61 gives 
that there exists such that G L^. The result now follows easily. □ 
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5. Studying (jCRp in the basis of special Hermite functions 

We are able to describe quite precisely the non-linear operator T, in the Hilbertian basis of L^(M^) 
provided by the special Hermite functions, which will be defined below. 

5.1. The special Hermite functions. We hrst recall some elements of the spectral theory of H 
and L. We follow mainly [U Appendix D], see also [23]. Dehne the creation and annihilation operators 


^{x + d^), al = ^{x-d:c) and ay'^= + dy), al = ^{y-dy). 


In the sequel, it will be more convenient to work in complex coordinates. Therefore we set z = x + iy, 
'z = X — iy and dz = ^{dx — idy), dz = \{dx + idy) before dehning 

CLd — {Ox iP'X ^^y') — 2 

We record the following formulas 

aS = ^(z-25^), al = ]^{z-2dz) 

a-d] ~ \P‘gi^g\ ~ ~ [o-gjad] = 0 

F{a*dU) = —ia*dU, F{a*u) = —ia*u 


(5.1) 


H — —Adzdz + 1^1^ — 2 ( 0 ^ ad + a*ag + l) 
L = zdz- zdz = a*d ad- a* ag 
X • V = zdz + zdz- 

We are now able to dehne the so-called special Hermite functions 

1 


and if n -|- m is even we set 
It is easy to show that 


'4’n,m — , -j-- {cid) (®g) 

^/^Tn\ ml 

n+m 71 —-m 


^ 1 




l,m? 

l'0n+l,ri 


which implies 
(5.2) 




n — m 


"V^rx—l,m+l 1 


^dVn,m — 


n + m 




(5.3) 




n — m -\-2 






n + m + 2 


"^n+l,m+l • 


It turns out that the families {il^n,m,n > 0,m > 0) and {ipn,m,n > 0,m > 0) are well-adapted 
to H, L and F. 
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Proposition 5.1. (i) The family {ipn,m-,n > 0,m > 0) is an Lf-normalised Hilbertian basis of the 

space L^(M^) such that 

Tdlfn,m — 2(71 -|- 7 Tl -|- — (n — 6 ^ ^ ^ ’f’rijm- 

(a) The family {ipn,m,n > 0,—n < m < n,nm even) is an L'^-normalised Hilbertian basis of 
such that 

TI^n,m — 2(71 + Lipyi Yti — '^Tn,rm Tn,m — 6 ^n,m- 

Proof. It follows in a straightforward way from the formulas (15.1|) . □ 

In other words, {ipn,m) and {ipn,m) are Hilbertian bases of common eigenfunctions of H, L and T, 
which is consistent with the fact that these linear operators commute with one another. 

Define the eigenspace En = span{ipn,Tn, —n < m < n,n -\- m even}. Let u G L^(M^) which can be 
written 

+0O n 

(5.4) U — ^ ^ rin; with Un — ^ ^ Cn^mTn,m G Eji, 

n=0 m=—n 

with the convention Cn,m = 0 if 77 + rn is odd. 

5.2. The operator T in the basis of special Hermite functions. 


m^dx ) 1/ ni+n2=n3+n4 I j 
I mi+m2=r7i3+m4 J 


,m4 


Proposition 5.2. We have 

(5.5) T{}Pn\,m\ iTn2,ra2iTn^,miiTnA,m4) — T I / Tni,miTn2,m2Tn3,m3Tn4 

\Jr^ 

and 

(5.6) T {}Pn\,mil Tn2,1712 1 ^11:^,1113) — ^ [I Tni,miTn2,m2Tn'j,m^Tn4,m4dx ) 

\Jr^ J 

with 774 = 77i + 772 - and 7774 = 7771 + 7772 - "73 (actually, / ipni,miTn2,m2¥i^ffmTTIfTfdx = 0 if 

J1R2 

7771 + 7772 m3 + m^). 

Proof. We first check that 

^{Tni,mi 1 Tn2,m2 1 Tnz,mz ) = £{}PrLi,mi 1 Tn2,m2 1 Tn^^mz 1 Tn^,1114)^114,,m^ 

with 774 = 77i + 772 — 773 and 7774 = mi + m2 — m3. Since (pn,m is a common eigenfunction of H 
and L with eigenvalues respectively n and 777, Corollary 12.51 implies that T{(pni,miiTn2,m2iTn3,m3) is 
an eigenfunction of H and L with eigenvalues respectively 774 and 7774. Thus it is collinear to ipn4,m4- 
The definition of £ gives the desired result. 

Next by (12.3p . 

TT 

^{.Tni,mi 1 Tn2,1112 1 Tn3,m3 1 T114,1114) — ® ^ ^ ^ dt j Pni,miTn2,m2Tn3,m3Tn4,m4dx, 

J-f yR2 

therefore if 771 + 772 = 773 + 774 and 7771 + m2 = m3 + 7774 

^i}Pni,mii Tn2,m2i Tn3,m3i Tn4,m4f — T / Tni,miTn2,m2Tn3,m3Tn4,m4dx, 

JM2 
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which concludes the proof. 

As a result, expanding / as in ()5.4p . the equation ()CE| is equivalent to 

iUn — ^ ^ T ^ Uyi2 1 , 


□ 


ni,n2,n3>0 

ni-\-n2—nz=n 


or: for all n > 0 and —m < m < n 


(5.7) = (f 

ni,n2,n3>0 —nj<mj<nj ^ ® t' 

711+712—n3=n 7711+7712—7713=771 

5.3. Conservation laws in the basis of special Hermite functions. 

Lemma 5.3. Let Q be one of the operators ad, ag, or a*. Then for all /i, / 2 , /s £ 5(]R^) 

g(r(/i, /2, /s)) = r(Q/i, /2, /s) + nh, q/2, /s) - n/i, /2, qVs). 

a consequence, 


ns,ms- 


/ (adw)u, / 

JR2 JK2 

are eonservation laws for dcEl). 


I, / (a»u, / (a>)a 


Proposition 5.4. ^et u = '^n=o'^rn=-n with the convention that Cn,m = 0 ifm + n is odd. 

In these coordinates, the conservation laws are the following 

+00 n 


-'71, m \ 


(5.8) 

(5.9) 

(5.10) [ 

Jr- 


/ (odw)! 

JR2 


ii^iil2(r2) - 

n=0 m=—n 
+CX) n 

II'^'^IIl2(r 2) = 2 ^^(n + 1) |cn^7 

77=0 m=—n 

„ +00 77 

(Ln)u = Yi E 

n=0m^-- 

^ +CX) n 

E E 


^ Cn.m 


/ ( 

Jr2 


77=0 m=—n 

+00 77 


n + m + 2 


<^77 + 1 , 777+1 Crfim 


anUju = / ua*u = 


I |2| |2 

z m = 


E E 

n=0 m=—7i 

+CXD 71 

^(n+1) \Cn,m\^ + Y, ^2 

71=0 rn=—7i n=0m=—n 


n — m + 2 


"Cn+l,m—1 ^ 71,771 


+ CXD 


\JrP‘ — 


(cn,m 2,m H“ ^n,7n Cn—2,m) 


+CX) n 


+CX) rx 


(5.11) f {x-Vu)u = - Y^ Y |c77,777p + X] ^ 


\Jn^ — 




■n,7n\ , 2 {pn,771 2,771 ^7i ,tti ^ti— 2,m) ' 

rx=0m=—n rx=0m=—n 


2 
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Proof. The three first relations are straightforward. For (15.811 and (|5.9I1 we use (15.21) and (15.311 . 
For (I5.10p , we write z = Ug + and z = Od + a* to get 


n — m 


^Pn,m — \ „ Pn—l,m+l T 


n + m + 2 


~Pn+l,m,+l 


and 


y/rP — m? y^(n + 2)^ — mf 


1^1 Pn,m — (ll' T 1)9^71,m T „ Pn—2,m T 


V^n+2,m • 


We now turn to ()5.11ll . Similarly, 


dzPn,m — „ 


1 n + m 


1 n — m + 2 


2V 2 2 V 2 


"V^n+l,m —1 ? 


and 


zdz^Pn,m = ^{m- - m?ipn-2,m “ ^\/(n + 2)2 - nPipn+2,m, 


and using that ^pn,m = Pn-m we get 

zdzPn,m = - l)¥5n,m + jV - m‘^<fn-2,Tn “ ^ \/(n + 2)2 - m?ipn+2,m, 


which gives the desired result. 


□ 


6. Dynamics on the eigenspaces of H 

Recall that, for N G N, Eyj is the A^-th eigenspace of —A + |xp, associated to the eigenvalue 
2N + 2. It is spanned by the eigenfunctions ^N,m, m G In = {—H, —N + 2,..., A — 2, N}, thus 
natural coordinates on En are provided by the {cN,m) (which we simply denote (cm) when the context 
is clear): if u G En, 

U — ^ ^ CyriPN,m- 

There are only two independent conserved quantities for the restriction of (ICRp to En- 
M{u) = ||u||^2 = |cmp and P{u) = J Lu-u= m|cmP, 

the associated symmetries being of course phase and space rotation: u i-A and u i—>• R\u. 

6.1. Dynamics on Eq. The eigenspace Eq is generated by the Gaussian ipQfi{x) = ■ For 

data uit = 0) = coy^o.O) the solution u{t) = c(t)(/?o,o is given by c{t) = e“* 2 hol 
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6.2. Dynamics on Ei. Write u = c\Lpi^i + By Lemma 12.21 

8{u) = TT^ f 

JR2 

Using that f = f = l/(47r), it is easy to see that the Hamiltonian reduces on Ei to 

8(u) = TT^ / + C_lV9l _i|^ 

2K2 

= ^(|ci|^ + 1c_i|^ + 4|ci|2|c_i|2). 

The equation (ICRh can be written idtc = ^|§ or in other words 

rzci = ^(lcip + 2|c_ip)ci, 

\ = ^(2|cip + |c_ip)c_i. 

It is now easy to integrate this equation: if uq = 

u(t) = c? exp ( - ^(|c?|2 + 2|c° + c° 1 exp ( - ^(2|c?p + |c° ip))v9i _i, 

thus in particular, every solution is quasi-periodic; even more, every solution is an M + aP wave, for 
some a. 

Even if the variational structure on Ei is fairly simple, let us record it before moving on to the 
more complicated situation on E 2 . Maximizers of 8 for fixed mass are the {cj} such that |ci| = |c_i|; 
and minimizers of 8 for fixed M satisfy ci = 0 or c_i = 0. These give rise to M-stationary waves. All 
the other solutions are non-trivial M + aP waves, which are maximizers of 8 for M + aP fixed, with 
|a|<|. Ifl 7 ^|Q;|>^, the maximizers degenerate and are given by ci = 0 or c_i = 0, whereas if 
OL = 1,-1 the constraint M + aP degenerates and the maximizers can be infinite. 

Thus, all the solutions can be obtained as extremizers, and are orbitally stable, in the sense that 
the moduli |ci(t)|, |c_i(t)| are stable with respect to perturbations of the data, uniformly in time, but 
angles are not. However, this can be seen directly, without resorting to variational considerations! 

6.3. Dynamics on E 2 . We have a good picture of the dynamics in the eigenspace £" 2 . 


6.3.1. Writing down the equation. Decomposing u = C2(p2,2 + co¥? 2 ,o + C-2^2,-2, we find as above that 
8{u) = TT^ [ |u|^ 

iK2 


TT 

4 


^1C2|'^ + ^|C-2|^ + |co|^ + 3|c2p|c_2p + 2|c_2p|co|^ + 2|c2p|coP + C2C_2Co^ + C_2C2c| 


while the conserved quantities of (jCRp read 


M{u) = |c_2p + |coP + |C2 
P{u) = |c2p - |c_2p. 
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One can check that the functionals £, M, P are in involution on E2, which has dimension 6; this 
makes the system completely integrable. Using the formulation it = , the equation reads 

ic2 = ^ [3|c2|^C2 + 6|c_ 2|^C2 + 4|co|^C2 + 

< icL 2 = ^ [ 3 |c_ 2 pc _2 + 6 |c 2 |^c _2 + 4 |copc _2 + 2qco] 
ico = ^ [4|copCo + 4 |c 2 pCo + 4 lc_ 2 |^Co + 4 c 2 C_ 2 Cij] . 


Rescaling time by r 

( 6 . 1 ) 


lOvrt and switching to the unknown function dj(T) = e this becomes 

id2 = —|d2p<^2 + 2|(i_2|^d2 + 2(i_2(io 
< id-2 = ~\d-2\^d-2 + 2\d2'\^d-2 + 2d2d^ 

idf) = 4d2d-2do. 


6 . 3 . 2 . The M- and M+aP- stationary waves. A computation gives all the M- and M+aP-stationary 
waves: 


0 , 0)6*/^* 

with z G C, 

and // = 

|Z|2. 


), z, 0)e*^* 

with z G C, 

and [1 = 

|Z|2. 


), 0 , z) with z G C. 




;, P, 0)e*^* 

with z, z' G 

; C, |z| = 

\z'\, and fj, = —|zp. 




[7 jP^+P 2 

/ —P^ 2 

\J 9 ^ 

^ with A,/ 3 i,/32 G M, and g, = 

- 8 \2 
- gA . 

(#"■. 



1 gi/it / 3 i, /32 G M, and // = 

-fA2 


(g) {d2,d-2,do) = R-acotiz, z', 0 )e^‘^^ with |2;| / \z'\ non-zero, w = 7 = j^L, [^,L , and a = 


' 2u) 


(h) {d 2 ,d- 2 ,do) = 


bl^+bT 

with x,y,2,/3i,/32,Ai,P G M, /r e 




j/2 ^ (8^2 _ 3^)2 ( _ ^ 

and finally e = ±1 if ^ < 0, e = =bi if ^ > 0. 

Before discussing the stability of these solutions, let us show how the system can be integrated. 


6 . 3 . 3 . Integrating the equation. Setting (recall that P = |(i2p ~ \d-2\ 
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the equation (| 6 . 1 I 1 becomes 


iD2 = |-D_2| D 2 + 2Z)_2-Do 
iD-2 = \D2\‘^D-2 + 2 Z) 2 -Dq 
IDq = AD2D-2D0. 

Setting now D 2 = Z )_2 = Bq = C'e*'^, the conservation of mass and momentum give the 

relations 


A = 


M + P-C^ 


B = 


M-P-C^ 


which allow to eliminate A and B and obtain the new equation 

C = AABC sin(a + /3 — 27 ) 


a = —B^ — 2 -^^ cos(a + /3 — 27 ) 

/3 = — A ^ — 2 ^^ cos(q: + /? — 27) 

7 = —AAB cos(a + /3 — 27 ). 

Switching to the new unknown ^ = a + (5 — 27 , we obtain the two-dimensional ODE 

C = d^BCsin^ 

i = - A ^ - B ^ -2 + ^ - cos^. 

As a two dimensional ODE, it can be fully understood (by plotting the phase portrait) and then we 
can deduce the behavior of the full system. We do not pursue this direction here. 

6.3.4. Orbital stability of the M- and M + aP-stationary waves. It can be deduced from variational 
considerations, as well as the reduction to a two-dimensional ODE that was just presented. We 
examine one by one the waves presented above. 

(a) is orbitally stable since it maximizes the angular momentum for fixed mass. 

(b) is orbitally stable since it minimizes the angular momentum for hxed mass. 

(c) and (d) are not stable, and there are actually orbits joining arbitrarily small neighborhood of 
the former and the latter. This can be most easily seen by considering the reduced system in 

C = 20 (l-C 2 )sinC 


(C, in the case P = 0. Assume that M = 1 , then 
waves under consideration correspond to the orbits 


^ = (C^-l)+4(l-2C^)cos.e 


and 


C = 1 

^ = —A cos ^ 


C = 0 

^ = — 1 -|- 4 cos ^ 

and an analysis of the phase diagram gives the desired conclusion. 

(e) is orbitally stable since it minimizes the Hamiltonian for hxed mass, as a lengthy computation 
shows. 

(f) is orbitally stable since it maximizes the Hamiltonian for hxed mass, as a lengthy computation 
shows. 

(h) is orbitally stable since it maximizes (for /r < 0) or minimizes (for /x > 0) the Hamiltonian for 
M -|- aP hxed - for a properly chosen a. 
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6.4. Dynamics on E^. Some M-stationary waves on Ej^i are given by all the with k ^ 1^. 

The waves (pN,N and (p^-N are orbitally stable as extremizers of P for M fixed. By analogy with the 
case = 2, it seems natural to expect that the other (pM,k are unstable. 

These are by no means the only stationary waves: for instance, it is easy to check that + 

z'pN-k gives rise to a stationary wave if \z\ = \z'\, and k,k' > y. Still by analogy with the case 
N = 2, this wave should also be unstable. 

Other orbitally stable waves should be obtained by minimization, or maximization of £ for M, or 
M + aP fixed. It seems plausible that this extremization procedure should produce new waves than 
the ones which have already been described. Focusing on the case where the mass is hxed: M = 1, 
this would be the case if 

(6.2) max £{(p) > max£{ipN k), respectively min £{(p) < min £{(pN k)- 

M{ip)=l fcG/jv ' M{ip)=l kGl]\^ 

This does not follow from known estimates: we only know that 

max £{ip)KiN~^, £{ipN,o) < N~^{logN), and £{ipn,n) = £{pn-n) ^ 

M{ip)=l 

(the first estimate is taken from [T7], the second one, valid for N even only, from m, and the last 
one is a simple computation). However, there are good reasons to believe that, for instance, the first 
inequality in ()6.2I) should hold. Indeed, the near maximizers of £{p), as explained in [T7j, are expected 
to focus along rays, or at points. This is not possible for the pN,k, which satisfy \(pN,kiz)\ = pN,k{\z:\)- 


7. Dynamics on the eigenspaces of L 
Adopting radial coordinates (r, 0), let us set 

Fn = {e^"^f{r), with/GL2(^2)|_ 

It is the n-th eigenspace of the rotation operator L, which is left invariant by the dynamics of (jCRp . 
An eminent instance is the set of radial functions Eq. 


7.1. Dynamics on Eq, the set of radial functions. The following can be found in [231 Chapter 1] 
and [23l Corollary 3.4.1]: the Laguerre polynomial L® of type 0 and degree k > 0 is defined by 


(7.1) 


^ k\ dx^ 


(e a; G M. 


These polynomials are orthonormal on T^([0,+oo),e ^dx) 

(7.2) j Lf\x)Lf\x)e~^dx = Sjk, 

and are related to special Hermite functions of second index 0 by 


^2k,o{^) = ^Lf\\x\‘^)e 
\/vr 


( 0 ), 


Simply denote hj^ = (p 2 k,o^ then Hh^ = (4A: + 

^(^ni ; ^77.2 ? ^77-3 ? ^71-4 ) — ^ ^ ^771^772^713^714)1771+712=713+774 




and 
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T{hni 1 hn2 1 hn:^) — TT ( ^n2^ 714 )^n4) ^4 — n-i -|-77-2 713 . 


+ 00 

Write / = Cn/in- Then by (15.71) the equation (jCRp is equivalent to 

n=0 


iCyi — TT ^ ^ Ctt,^ C7^2 

ni,n2,n3>0 

It was already established in [7] that centered Gaussians generate stationary waves in Fq. Another 
stationary solution exhibited there is the self-similar function which does not belong to but is 
in the generalized 0-eigenspace of L. Finally, all the hn give rise to stationary waves: u(t) = e~^^”*hn, 
with OJn = S{hn). 

Other stationary waves can be obtained by letting the symmetries of the system act on them. The 
symmetries of dell) which leave the set of radial functions invariant are u i—)• u i—)• 

u 1 —^ u f—^ S\u = Am(A-) and u i—)■ Since the lens transform formula ()2.5p expresses 

the fifth symmetry in terms of the four first symmetries, it suffices to consider the four first ones. 
Applying them to hn gives the orbit 

On = ( 0 , n, fi, A) e M^}. 



In the case of the Gaussian, we obtain the orbit 

00 = {0, O h, A) G M^}. 

It was proved in [7] that the Gaussian is orbitally stable in in the sense that data close to 

in yield solutions remaining close to Oq for all later times. 

For this reason, it is interesting to express the orbit Oq in the (c„) coordinates: 

, hn) = , K). 


We now successively change variables to z = |a:p, use the formula (17.ip giving hn, and integrate by 
parts repetitively to obtain 


(e 




SxhQ , hn) — 

TT 






A / d 

n! Jo \dz 

(-l)’^A /I-A2 


n\ 


+ 


^+00 

i ' 




1-A^ 


-I- 7/7 


-hi/i 


A 


n+1 * 
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Therefore, S\hQ , hn) 

coordinates 


_pi0 iv{in+2) \ V 2 

^ ^ ^ ( 14-X2 , . 


and the orbit Oq reads, in the (c„) 



g*6g*i'(4n+2)^- 


+ i^J. 


1+A2 


+ inj 


\ n+1 




7.2. Dynamics on Fn, n 7 ^ 0. The subspace does not contain Gaussians anymore, and the role 
of the ground state is played by 





It is not clear whether it minimizes £ on Fn for fixed mass; but it is clearly the minimizer of {Eip, cp) 
on Fn for fixed mass. This gives immediately orbital stability on Fn for the M-stationary wave arising 
from (pn,n- We study the dynamics on this family {p}n,n} in the next section. 


7.3. Other eqnivariant stationary waves. We first claim that the equality 




^inO 

= OJ - 

r 


holds for some real number uj (recall that T is bounded on such that rf G thus the 

left-hand side makes perfect sense). Recall that F commutes with the rotation operator L and the 
dilation operator S\. Therefore, must be invariant by S\, and of the form e*"'®/(r); 

thus it has to be equal to for some co. 

As a consequence, we obtain a new stationary wave in the generalized n-th eigenspace of L: —. 


8. Dynamics on the Bargmann-Fock space 

Denote by 0(C) the space of the entire functions in the complex plane. Then the Bargmann-Fock 
space is given by n ( 0 (C)e“hl /2^_ admits an orthonormal basis given by the special Hermite 

functions (pn,n, which we will simply denote p>n to alleviate notations: 

(8.1) ipn{xi,X2) = ^L=(xi -h 

V vrn! 

Recall that ipn is such that Hipn = 2(n -|- l)(pn and = 1. 

It is an invariant subspace for (ICE,p . and we consider in this section its dynamics restricted to it. 
Out of all the symmetries of (ICRD . only three act on its restriction to the Bargmann-Fock space: 
phase rotation u e->■ rte*®, with 0 G M; space rotation u e-)■ Rqu, with 0 G M; and magnetic translations 
u eA + for ^ G C. These symmetries are associated by Noether’s theorem to the three 

conserved quantities 

M= f \u{z)\‘^ dz, P = [ (kP — l)| 4 t(- 2 )P Q= f z\u{z)\‘^dz. 

Jc Jc Jc 
















27 


8.1. The framework. 

Lemma 8.1. Let ipn be defined by (|8.ip . Then 

T (v^ni 1 Tn2 ) ) — Q^ni ,n2,n^,rnTni j 714 — '^1 T 712 713, 

with 


( 8 . 2 ) 


a. 


vr (771 + 712)! 


ni,n2,n3,n4 


— £{Tnn Tn2i Tnsj Trii) — /-j-j-j-r Ini+n2=n3+n4 • 

2 2^i+”’2.^4j^!7j2!r73!n4! 


Proof. The first claim follows from Proposition 15.21 

With the change of coordinates 2 ; = re*® and p = 2r^, we get 


£{(p n\i Tn2i Tn^i Tnfij — ^ / +ni+n2+113+^4 

2 R 2 

= , . \ / (xi+7X2)'^l+"^(xi-7X2)’^^+"^e-2|-l'dx 

Vni 17721773 1774 ! Jm 2 


r+co 


^ _ f g7(ni+n2-n3-n4)0^^ / ^ni+n2+n3+n4+lg-2r2^^ 

1772 !7731774! Jo 


vr 


+ 00 


2.2”'l+”2.^ni!772!773!774! V 7o 
7 r( 77 l + 772)! 


^ni+n2e- 


ni+n2=n3+n4, 


2.2”i+”2y7j.j^]j.j2!r73!774! **^^**^ «3+>T’4; 

which was the second claim. 

As a consequence we a have the following result. 

Lemma 8.2. Denote by 11 the orthogonal projector on the space . Then 

/ X 1 |z|^ f — + 1 ^ 

(nu)( 2 ;) = —e 2 J 2 u{w)dw, 

and we have 

T{iPni,Tn2,Tn3) = TT‘^D.{TniTn2¥^) ^+ni (6+n2 

Proof. The hrst point follows from the fact that the kernel AT of 11 is given by 


□ 


=-ldV 2 


+00 .. 

^—•' vr 

n=0 

The reformulation of T in terms of 11 is then a direct computation using a polar change of variables. □ 
Hence the (jCR| equation reads on L^(]R^) n (C>(C)e“l^l^/^) 

idtu = 7rn(|77p77). 

We remark the resemblance of this equation the Szego equation of Gerard and Grellier where there H 
is the Szego projector [9]. 
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8.2. Stability in of stationary waves. Consider a Gaussian solitary wave in the Bargmann-Fock 

up 

space: u{z) = -^e ~Its stability follows directly from the results in [7]. 

• It is orbitally stable with respect to perturbations in (ggg |7J Proposition 6.8]). 

• With respect to perturbations in L^, it is orbitally stable modulo the symmetries acting on 

the system (phase rotation, space rotation, and magnetic translations). In other words, any 
perturbation of it remains close to with 0 G M and ^ G C} (see [71 

Proposition 8.5]). 

We now would investigate the stability of the stationary waves with un = £{(Pn)-, for 

> 1. Since nonlinear stability seems to be a delicate question, we focus on linear stability (under¬ 
stood as the absence of exponentially growing mode). 

Proposition 8.3. For N >{), consider the wave . The following stability/instability results 

hold in the Bargmann-Fock space: 

(i) If N = 0 or N = 1, this wave is linearly stable. 

(a) If N >2, it is linearly unstable. 

(in) The number of unstable modes of the N-th wave is o{N). 

Proof. We already have the result for the Gaussian, thus we can assume that > 1. It will be 
convenient to use the basis provided by the special Hermite functions ipn- write u = Yl^=o^nTn- The 
stationary waves whose linear stability we will investigate read in these coordinates 

Cn — 

(where 6 is the Kronecker delta function). The linearization of (ICRp restricted to the Bargmann-Fock 
space around this stationary wave is given by 

(8.3) icN = 2aNNNNCN + 

(8.4) ick = 2akNkNCk + C2N-k if A: G {0 ... A^ — 1} U {N -|- 1... 2A^} 

(8.5) ick = 2akNkNCk k >2N -\-l, 

where a is given by (18.2p . The equation (18.51) is obviously stable. As for the equation (18.31) . the change 
of unknown variable cv = leads to the equation 

iz = aNNNN{z + z), 

whose solutions grow at most linearly. Finally, to study the equation (|8.4p . observe that it only couples 
Ck and C2N-k- Set A: G {0 .. . A^ — 1} U {A^ -|- 1... 2A^} and write 

U = OnNNN, Ck = C 2 N-k = C~'^‘^^y, A = ak, 2 N-k,N,N, B = OkNkN, C = a 2 N-k,N, 2 N-k,N- 

The equation (18.4p becomes 

( ix = Ay -|- {2B — uj)x 
\^iy = Ax + {2C - u})y, 

which implies 

X i{2B — 2C)x — (A^ — {2C — u}){2B — uj))x = 0. 
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This equation has exponentially growing modes if and only if its discriminant is positive: 

A{N, 2N-k) = A{N, k) = 4(^2 - {B + C - ujf) > 0, 
for 0 < A: < — 1. The results of the proposition are then implied by the following lemma. □ 

Lemma 8.4. With A{N,k) given by the above definition, 

(i) A(1,0) < 0. 

(ii) A{N, N -2)>0 for all N >2. 

(Hi) For any A € (0,1), A{N, /c) < 0 for k = XN and N sufficiently big. 

Proof. The part (i) is straightforward. We then check (ii). Thanks to the formula giving a, we find 
iov k = N — 2 


A + B + C — LO = aN-2,N+2,N,N + 0:n-2,N,N-2,N + 0'n+2,N,N+2,N “ 0:n,N,N,N 

'k{2N-2)\ / 2(2iV-l) (4iV2 _2jV-5)(A + l) \ 

“ {N - 1)\{N - 2)!22Af+i \ y(]v + 2)(iV + l)A(iV^ + (iV + 2)(iV + l)iV(iV - 1) J ' 

Since 4A^^ — 2N — 5 > 0 for A > 2, we find that A + B + C — uj > Q ioi aA N >2 Next, 


A — B — C + UJ = aN-2,N+2,N,N — aN-2,N,N-2,N “ OIn+2,N,N+2,N “ 0']S[,N,N,N 

7r(2iV-2)! / 2(2iV-l) {AN‘^ - 2N - 5){N +1) \ 

~ (N - 1)\{N - 2)!22^+i Vyp + 2)(iV + l)iV(iV - 1) " (iV + 2)(iV + l)iV(iV - 1)J ' 

This term is positive iff 

4(2iV - l)2(iV + 2)N{N - 1) > (4A^^ - 2N - 5)'^{N + 1), 

or equivalently 8N^ + 52A^ — 53N — 25 > 0, which is the case for all N > 2. Since A = 4:{A + B + 
C + lv){A — B — C — cu), we find that A{N, N — 2) > 0 for all N > 2. 

To prove (Hi), we use Lemma 18.11 and the Stirling formula to prove that, if A: = XN for some 
A G (0,1), as N —^ oo. 


/ vr 1 

( ^ 

/ v^A(2 -X)VN 

^At(2-A)i-t^ 


B 


C 


UJ 


X 


3-A 1 


Observe now that, if A G (0,1), 


-/tt 1 
1 

A^(2-A)^-^ 


l + A 1 /(1 + A)i+^ 

Tiv V 


N 


AA21+A 

(3 - A)3-^ 


2-Xy/N V(2-A)2-^23-a 


N 


N 


(l+A) 


l+A 


’ ^’^d also belong to (0,1). This 

implies that A, B, and C decay exponentially in N while uj only decays polynomially in N. Since 
A = 4(^2 — (^B + C + n;)2), this implies in turn that, for A fixed, A{N, XN) < 0 ioi N sufficiently 
large. □ 


























30 


PIERRE GERMAIN, ZAHER HANI, AND LAURENT THOMANN 


8.3. The framework. It is easy to check that the orthogonal projector on 0{C)e is bounded 
on L°°, and thus can be naturally extended from to Therefore, the equation 

idfU = 7rn(|rtpM) 


is locally well posed in for data in L°°. 

Here we make the link with a result of Aftalion-Blanc-Nier [H Theorem 1.4]. Denote by Hhoi the 
orthogonal projection on the space of entire functions on C (in [T] this corresponds to H/^ with h = 1). 
For /, set u = then 


nu = n(e-i^i'/2/) 





Thus 

n(|n|2u) =e-l^l'/XoKe"'^'VlV). 


As a consequence, the function Uriz) = e given by [U Theorem 1.4] is a stationary solution 

to ()CB| . The function fr belongs to the space 


Appendix A. Two general structures 

The various equations derived in the present paper present striking similarities, they all belong to 
one, or two, of the general structures described below. The presentation we give is only formal. 


A.l. An equation on sequences. Assume that (o-kimn) (k,i,m,n)eA'^ ^ where A C Z satisfies the sym¬ 
metries 

O^klmn — Oilkmn and OLklmn — Olmnkl- 


Then the equation 


iCn 


E 


Q^nlr^ 27 ^ 3 r^Cnl ^722^^3 


721 + 722 = 723+71 


derives from the Hamiltonian 


GiiCn)) 



721+722=723 + 724 


^ 721722723724^721 


Conserved quantities for this equation are the norm as well as the Hamiltonian. 

Examples are 

• A = Z, anin 2 n 3 n 4 = Ij which is simply the equation iii = \u\‘^u on T seen in Fourier space. 

• A = N, anin 2 n 3 n 4 = lfc,«,m,n>o, which is the Szego equation. 

• A — {0,..., A}, cxnin 2 n 3 ni — £‘pN,n 2 iii which is the equation on En- 

• A = N, Q:nin 2 n 3 n 4 = f 2 ”i+” 2 Vm?n 2 !n 3 !n 4 ! is the Lowest-Landau-Level equation. 

• A = N, anin 2 n 3 n 4 = + / ^^ 2 + 3 where the (hn) are the normalized radial Hermite 

functions, which is the ’’radial equation”. 
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A. 2. An equation on functions of a continuous variable. Assume that E is a closed subset of 
with orthogonal complement E-^. Let If be the orthogonal projection on E. 

The equation (whose dependent variable is a function u on M^). 

iu = Yi {\u\^u) 


derives from the Hamiltonian 

£{u) 


Conserved quantities are / |u|^ and / |up. 
Instances are 


u\^ on E. 


• E = which is simply the equation ii = 

• k = 1, E = (functions with positive frequencies), which is the Szego equation. 

• k = 2, E = E]\[, which is the equation on E]\[. 

• k = 2, E = /^0(C), which is the Lowest-Landau-Level equation. 


Acknowledgement: The authors are grateful to Patrick Gerard for pointing out the connection 
with the Lowest-Landau-Level equation. This work was initiated during the visit of the third author 
to the Courant Institute of Mathematical Sciences, and he thanks the Institute for its hospitality. 
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